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ABSTRACT 

Sub-arcsecond scale Keck images of the young A1V star, 49 Ceti, resolve 
emission at A = 12.5 and 17.9 /im from a disk with long axis at PA 125 ± 10°and 
inclination = 60 ± 15°. At 17.9 /im, the emission is brighter and more extended 
toward the NW than the SE. Modeling of the mid-infrared images combined 
with flux densities from the literature indicate that the bulk of the mid-infrared 

t— i \ emission comes from very small grains (a ~ 0.1 //m) confined between 30 and 60 

AU from the star. This population of dust grains contributes negligibly to the 
significant excess observed in the spectral energy distribution. Most of the non- 

\/~} [ photospheric energy is radiated at longer wavelengths by an outer disk of larger 

grains (a ~ 15 /im), inner radius ~ 60 AU, and outer radius ~ 900 AU. Global 

O \ properties of the 49 Cet disk show more affinity with the (3 Pic and HR 4796A 

disks than with other debris disks. This may be because they are all very young 
(t < 20 Myr), adding strength to the argument that they are transitional objects 
CLf between Herbig Ae and "Vega-like" A stars with more tenuous circumstellar 

O ' disks. 

a: 

Subject headings: circumstellar matter - infrared: stars - planetary systems: 
.^j | formation - planetary systems: protoplanetary disks - solar system: formation 

- star: individual(49 Cet) 

1. Introduction 

Vega-excess stars were identified in observations with the Infrared Astronomical Satel- 
lite (IRAS) as having infrared excess emission attributed to circumstellar dust (see reviews 
by Backman & Paresce 1993 and Zuckerman 2001). More than a hundred Vega-excess candi- 
dates have been identified by comparison of IRAS sources with stellar counterparts in optical 
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catalogs since the early discovery of prototypes such as Vega and (3 Pic (Backman & Paresce 
1993, Mannings & Barlow 1998, Silverstone 2000). Estimated ages and circumstellar dust 
masses support the view that these objects are largely distinct from their pre-main sequence 
counterparts, the T Tauri and Herbig Ae/Be stars. The latter have typical ages less than 10 7 
yrs and fractional infrared luminosities, r = Li r /Lboi, of approximately 0.25 (Adams, Lada 
& Shu 1987), while Vega-excess stars may be as old as 10 9 yrs with values of r typically in 
the range 10 -4 to 10~ 7 . Larger values of r = 10~ 3 are found for a few Vega-excess stars with 
much younger ages (t ~ 10 7 yrs) that often appear to be associated with young kinematic 
moving groups (Moor et al. 2006). Jura et al. (1998) noted three A-type stars in the Bright 
Star Catalogue with r ~ 10~ 3 and interpreted them as representative of an early phase of 
debris disk evolution. f3 Pic, HR 4796A, and 49 Cet have since been assigned stellar ages of 
20 Myr or less (Stauffer et al. 1995, Zuckerman et al. 1995, Zuckerman et al. 2001) and may 
well be important as examples of the transition from viscous protostellar accretion disks to 
planetary systems. The existence of such a phase is supported by recent Spitzer studies of 
weakline T Tauri stars in nearby clouds where a small sample of objects with values of r in 
the high debris disk range was recently discovered by Cieza et al. (2006). 

High-resolution images can reveal features that signal the presence of gravitational in- 
teractions between unseen planets and orbiting dust grains. In recent years more than a 
dozen such disks have been resolved, revealing tantalizing features like ring arcs, clumps, 
disk warps, and other non-axisymmetric signatures of potentially planetary importance 
(e.g. Holland et al. 1998, Koerner et al. 1998, Schneider et al. 1999, Heap et al. 2000, Ko- 
erner et al. 2001, Wilner et al. 2002, Clampin et al. 2003, Wahhaj et al. 2003, Ardila et al. 
2004, Liu et al. 2004, Schneider, Silverstone & Hines 2005, Wahhaj et al. 2005, Kalas et 
al. 2005 and Kalas, Graham & Clampin 2006 ). Simple gaps and rings are often apparent, 
although these may not be definitive signatures of unseen planets (Takeuchi & Artymowicz 
2001, Klahr & Lin 2001). Some features, however, are as yet unaccounted for by any other 
means. Examples include stellar offset from the orbital center of an elliptical ring around 
Fomalhaut (Kalas et al. 2005). The evidence for a planet-disk connection is now sufficiently 
compelling that the Vega-excess phenomenon has been proposed as an indirect detection 
method for planetary systems (Zuckerman & Song 2004). 

Dynamically important disk features have been identified in high-resolution mid-infrared 
images for 2 out of the three young A-type stars mentioned by Jura et al. (1998). Keck studies 
of HR 4796A (Koerner et al. 1998; Telesco et al 2000; Wahhaj et al. 2005) and f3 Pic (Wahhaj 
et al. 2003; Weinberger et al. 2003; see also Gemini mid-infrared observations by Telesco et 
al. 2005) reveal gaps and rings with azimuthal asymmetries, inclination offsets, and radial 
changes in grain size and distribution. Here we present Keck mid-infrared imaging of 49 Cet, 
the third nearby A-type star identified by Jura et al. (1998) as having a young debris disk in 
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an early transitional phase. 49 Ceti is an A1V star with an age of 8±2 Myr (Zuckerman et al. 
1995; Thi et al. 2000) and a Hipparcos distance of 61±3 pc. As discussed by Jayawardhana 
et al. (2001), mid-infrared radiation of 49 Cet (0.38 Jy in the IRAS 25 /im band) is reduced in 
comparison to HR 4796A, posing an added challenge to imaging at mid-infrared wavelengths. 
Nevertheless, these authors report radiation from the circumstellar dust that is marginally 
resolved at 10.8 /im with flux density 250 ± 50 mJy Substantial long-wavelength excess has 
also been detected at A = 60 /im by Sadakane & Nishida (1986) and at A = 1.2 mm (-F 12mm 
= 12.7 mJy) by Bockelee-Morvan et al. (1994). Here we present Keck II images of 49 Cet at 
A = 12.5 and 17.9 /im together with simultaneous modeling of both the images and spectral 
energy distribution. 



2. Observations and Results 

49 Cet was imaged with the Mid-InfraRed Large-well Imager (MIRLIN; Ressler et al. 
1994) at the f/40 bent-Cassegrain focus of the Keck II telescope on UT dates 7 October 1998 
and 28-31 July 1999. At Keck, MIRLIN has a plate scale of 0".138 pixel" 1 and a 17". 5 field of 
view. Filters centered on A = 12.5 and 17.9 /mi with widths 1.16 and 2 /im were used for the 
observations presented here. Initial sky subtraction was carried out by differencing images 
chopped 7" in an east-west direction at a frequency of 4 Hz. The telescope was nodded the 
same distance in a north-south direction after co-adding several hundred frames at each of 
the chop positions. The total integration times were 20 and 66 minutes at 12.5 and 17.9 /im 
respectively. Residual background emission in the form of a striping pattern was removed by 
masking on-source emission and separately subtracting the median emission within each row 
and column of pixels. Disk emission in each of the double- differenced chop and nod frames 
was shifted and added to produce a final combined image. Observation and data reduction 
of standards a Boo, a Tau, x Get and Vega proceeded in the same way for photometric 
calibration and served as a representation of the Point Spread Function (PSF). Images of 
49 Cet and PSF stars from separate nights of observation were co-added after weighting by 
the signal-to-noise ratio of their peak flux. Seeing effects on low-level emission in the PSF 
stars were mostly evident at 12.5 /tm and were reduced for both the source and calibration 
images by smoothing with a circular hat function. 

Final processed images of 49 Cet at angular resolutions of 0".41 and 0".48 at A = 12.5 
and 17.9 /im, respectively are shown in Fig. [TJ The original 64 x 64 pixel (8.8" x 8.8") 
images were finely gridded to 256 x 256 pixels and smoothed by convolution with a circular 
hat function of diameter equal to the Full Width Half Maximum (FWHM) of the PSF star a 
Boo. Images of a Boo, resampled and smoothed in the same manner, are shown in the insets 
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in Fig. [TJ Images used later for model-fitting purposes were smoothed in the same way, but 
not resampled beforehand. In Fig.[T]both images of 49 Cet are elongated along approximately 
the same position angle. To characterize this effect more accurately, normalized strip cuts 
were obtained from the images by windowing the emission in rectangular apertures aligned 
with the long axis of emission (PA = 125°) and of width corresponding to the PSF FWHM. 
The results, displayed in Fig. [2J, demonstrate that 49 Cet's emission profile is clearly resolved 
at both wavelengths. 

Flux densities for 49 Cet were derived using aperture photometry of images of the 
calibrator stars. Calibrator flux densities were adopted from the ISO Calibration web site 
( |http: / / www.iso.vilspa.esa.es/users/expl_lib/ISO/wwwcal/[ ) as follows: a Boo (-F12.5 = 478 
Jy; F 17 . 9 = 232 Jy), a Tau (F 12 . 5 = 423 Jy; F 17 . 9 = 205 Jy), X Cet (F 12 . 5 = 0.68 Jy; 
Fn, 9 =0.33 Jy) and Vega (F 12 . 5 = 24.9 Jy; F17.9 = 12.2 Jy). The resulting flux densities 
for 49 Cet were 200±26 mJy and 186±25 mJy at 12.5 and 17.9 /an. Although our 12.5 
fim value and the appreciably higher IRAS 12 /im measurement (330 mJy) are at first sight 
inconsistent, the discrepancy is explained readily by the fact that the IRAS band is much 
broader (6 to 18 /an) and thus includes more of the bright photospheric emission at shorter 
wavelengths. Peak brightnesses, the flux in the brightest pixel divided by the pixel area, were 
231 mJy/arcsec 2 and 583 mJy/arcsec 2 at 12.5 and 17.9 /im respectively. The corresponding 
rms background values were 41 mJy/arcsec 2 and 26 mJy/arcsec 2 . The extent of the long 
axis enclosed by the 2a contours of the 12.5 /im emission is ~ 1".64 (100 AU). At 17.9 /xm, 
elongation at the 2a level is ~ 2". 3 (140 AU) with the SE side clearly more extended than 
the NW side. 

Available flux densities for 49 Cet at mid-infrared and longer wavelengths are listed in 
Table 1. A 200 /im flux density from ISOPHOT observations (Walker & Heinrichsen 2000) 
was suspiciously low (0.32 Jy) compared to the value implied by an interpolation between 
the 100 /xm IRAS and 1.2 mm IRAM observations. We have therefore listed only the values 
obtained directly from the Infrared Space Observatory (ISO) Data Archive (1.1 ± 0.5 Jy at 
170 /im and 0.75 ± 0.5 Jy at 150 /im). Photometry in U,B,V,J,H,K,L and M bands from 
Sylvester et al. (1996) were used to constrain stellar photospheric properties as outlined 
below. 



3. Modeling 

We use a simple model of a flat optically thin disk with relatively few parameters 
(Backman, Gillett & Witteborn 1992; Koerner et al. 1998; Wahhaj et al. 2003; 2005) to 
derive the properties of the 49 Cet disk, since our images have limited spatial resolution, and 
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there are relatively few flux density measurements. We assume that thermal radiation from 
an annulus of width dr at a distance r from a star is given by: 

ft(r) = a(r)s x B[T p (rlX](^), 

where a(r) is the fractional surface density, T p (r) the grain temperature, B[T p (r),\] the 
Planck function, and D the distance. For 49 Cet, D — 61 pc. Moderately absorbing 
dielectrics with an effective grain radius, a, have radiative efficiency s\ = 1.5a/ X for A > 1.5a 
and E\ = 1 for A < 1.5a (Greenberg 1979). The grain temperature for efficient absorbers and 
inefficient emitters is T p (r) = 432a~^ 2 (L*/L ) a2 (rAu) -0 ' 4 K can be calculated from radiative 
balance equations. For very small grains, a < 0.05 /mi, which are both inefficient absorbers 
and emitters, T p (r) = 636(L ! „/L ) 2 / 11 (rAu) _4/11 (T*/T ) 3 / 11 K (Backman & Paresce 1993), 
where L* is the stellar luminosity, 18.4 L for 49 Cet. The two temperature laws converge at 
a = 0.05 /iin, and varying the grain radius across this point does not result in a discontinuity 
in temperature. However, the latter law yields degenerate solutions between <r(r) and a, 
and so only the product of the two quantities can be constrained. This degeneracy finds 
natural expression in the probability distributions generated by the Bayesian method which 
we discuss later. The assumed power-law emissivities are admittedly crude approximations 
for the grain temperature behavior. However, the constraining power of the extant data 
is also limited. Thus, more realistic temperature laws are counter-productive since they 
unnecessarily complicate the modeling method. The modeling constraints yielded on grain 
radii suggests a range of temperatures. It is better to work from these temperature ranges 
and investigate what kinds of dust grains they are consistent with, under more realistic 
modeling assumptions. 

To describe the disk morphology and orientation we introduce parameters Ri n , the 
inner radius, AR , the width, 7 , assuming radial density profile, a ~ r~ 7 , 0, the inclination 
of the disk to line of sight (90° is edge-on) and 9, the PA of long axis. We also include an 
asymmetry parameter f3, the factor by which emission from the NW ansa exceeds that from 
the SE. We simulate the stellar photosphere by fitting a Kurucz (1993) model of a star with 
T e //=9250 K to the optical and infrared photometry of Sylvester et al. (1996), and use the 
mid-infrared and longer wavelengths estimates from Table 1. The resulting photospheric flux 
density is added to the central pixel of the model image, and the final image is convolved 
with the appropriate PSF star to construct the simulated emission map. Models are varied 
over all parameter space, and the probability of each model calculated from x 2 following 
Wahhaj et al. (2005). 

Initially, the two Keck/MIRLIN images were treated separately in the model- fitting 
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process to differentiate their separate roles in constraining disk parameter values. Orientation 
parameters like <fr and 9 can be estimated only by fitting to these individual images. Most of 
the emission appears to arise in a small inner disk of very small (a < 1 /im) hot grains that 
contribute neglibly to the spectral energy distribution (SED). Both the temperature and the 
size of grains are mainly constrained by the relative brightness of the disk at A = 12.5 and 
17.9 /im. 

Our subsequent efforts to fit the SED and mid-infrared images simultaneously with a 
single disk model failed to reproduce the extensions in emission in the mid-infrared images. 
Unreduced \ 2 values for the best fit were ~ 45, with cardinality of the data only 27 (see 
also in Wahhaj et al. 2005). Fits to the SED alone required a much larger disk than that 
mandated by the images to reproduce the long wavelength flux densities and yielded only 
weak constraints on R in and grain radius a. We then invoked a 2-part disk model comprising 
an inner disk of small grains to simulate the MIR images and a large outer disk with a 
population of large dust grains to fit the SED. To introduce as few parameters as possible 
the inner radius of the outer disk is set to begin at the outer edge of the inner disk. These 
two data sets are complimentary in that the 17.9 /im image constrains the inner radius of 
the outer disk, while the SED limits the inner radius of the inner disk. This model effectively 
reproduces all the data, and yields a \ 2 of 28. Most probable parameter values for the inner 
disk are (p = 55 ± 20°, 9 = 120 ± 15°, R in = 30 ± 10 AU, AR = 30 ± 20 AU, a < 1.0 /im, 
a = 0.1 — 1.5 xlO -3 , (3 = 0.5—2, with no constraints on 7. The probability distribution for the 
grain radius, a actually peaks around 0.1/im. The distribution suggests no lower limit mostly 
because the temperature law becomes degenerate between grain radius and optical depth for 
a < 0.5/iinas discussed earlier. Outer disk properties are = 60 ± 15°, 9 = 125 ± 10°, Ri n = 
60 ± 15 AU, AR = 900 ± 400 AU, a = 15 ± 10 /im, a = 6.5 ± 1.5 x 10" 4 , /3 = 1 - 10, and 
7 = ± 0.5. There exists an unrefereed publication of an 850/zm JCMT/SCUBA detection 
of 49 Cet (Song, Sandel & Friberg 2004). At 8.2±1.9mJy, this value is inconsistent with the 
1.2mm IRAM detection. We note that if the IRAM flux were ignored and instead the JCMT 
flux were used, the outer disk width would have to be ~ 300 AU. The weak constraints 
on the disk asymmetry parameters can be explained by the fact that the average signal 
from the outer regions of the disk is roughly 3cr. Thus any asymmetry seen in our images 
are statistically not very significant. The probability distributions for the inner/outer disk 
configurations are dispayed in Fig. [5] & [6] . 

Images simulated from the two-part disk model fit to Keck images and the SED for 
49 Cet are displayed in Fig. [3J The simulated SED is shown in Fig. HI We recognize that 
even this model may be an approximation to some more complicated structure that is not 
discernable from the current data set. Nevertheless, a number of useful properties emerge 
: 1) the 2-component disk has a uniform orientation, within uncertainties; 2) characteristic 
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temperatures derived from a grain size-temperature relation (see Backman & Paresce 1993) 
are ~ 360 K for inner dust grains and ~ 95 K for the outer; 3) the total cross-sectional area 
of grains in the inner disk is ~ 5 AU 2 , and in the outer disk, ~ 1650 AU 2 ; 4) a grain density 
of p = 2500 kg m~ 3 (Jura 1998), leads to masses of 6.5 x 10 -6 M ffi and 0.32 M e respectively. 

4. Discussion 

Mid-infrared images have also identified variations in radial structure in the f3 Pic and 
HR 4796A disks. Structures with sharply varying density and grain size are apparent within 
the depleted inner region of (3 Pic's disk (Lagage & Pantin 1994; Backman, Gillett & Wit- 
teborn 1992; Wahhaj et al. 2003; Weinberger et al. 2003; Okamoto et al. 2004; Telesco et 
al. 2005). Likewise, a comparison of mid-infrared and Hubble images of HR 4796A demon- 
strates that it is surrounded by a narrow ring of large grains embedded within a wider ring 
of smaller grains (Wahhaj et al. 2005), properties that are consistent with dispersal of an 
exo-Kuiper belt by radiation pressure. The analysis presented here suggests a similar mor- 
phology for 49 Cet. It is evident in Fig. [5] that there is only a very small probability that the 
disk around 49 Cet extends all the way to the star. Indeed, it is most likely that the region 
interior to 20 AU is strongly depleted of dust. The absense of silicate emission in Spitzer 
Space Telescope observations of 49 Cet suggests that the circumstellar dust must be cooler 
than 200 K (Kessler-silacci et al. 2005). In combination with the constraints on grain size 
obtained here, this further strengthens the case for a large inner clearing of the disk. 

The major contribution to the excess luminosity comes from a large-grained outer disk 
with an inner edge that can barely be discerned in the 17.9 /xm image. Both the 12.5 and 17.9 
fim images predominantly trace a region that is approximately 30 to 60 AU from the star and 
composed of very small dust grains, a ~ 0.1 /mi. How can dust grains as small as 0.1 /im 
survive in the harsh radiation environment of 49 Cet? A conceivable explanation can be found 
in a paper by Artymowicz (1997). Grain-grain collisions result in a size distribution, n ~ 
q-3.5 (Dohnanyi 1969). Grains smaller than ~ 7 /im, similar to the blow-out (/3-meteoroids) 
size for HR 4796A (Wahhaj et al. 2005), exit the system under radiation pressure. However, 
as they travel outward along the disk mid-plane, they undergo further collisions producing a 
dust avalanche of even smaller grains. The resulting amplification in the dust production rate 
depends on the disk surface density and the width of the disk traversed. The smallest grains 
(called 7-meteoroids) , will experience a relative pile-up, since their outflow is braked more 
efficiently by residual gas (Thebault & Augereau 2006). An approximate expression for dust 
production-amplification due to the avalanche phenonemon is e N<Tm , where N is the number 
of particles produced per collision and a m is the fractional surface density experienced by 
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a particle traveling along the disk midplane. Collision physics predicts values between 10 2 
and 10 3 for N. Assuming that the disk thickness is r/10, and the inner disk is 30 AU wide, 
a m = 8 x 10~ 3 . Thus the amplification factor ranges from 2 to 3000. Detailed numerical 
simulations of the avalanche mechanism show that, in a debris disk with properties similar to 
49 Cet, an amplification factor of 200 may be expected (Grigorieva, Artymowicz & Thebault 
2007). Obviously, too high a dust production rate predicts too short a disk lifetime. At the 
same time, rates consistent with the pile-up of 7-meteoroids are accomodated by the range 
of predictions. 

Another plausible mechanism for 7-meteoroid production is the sublimation of icy bodies 
analagous to comets or Kuiper Belt Objects in our own solar system. Submicron refractory 
particles may be held together by volatile ices to produce large particles unaffected by ra- 
diation pressure. Upon sublimation, however, the smaller particles may be released and 
subsequently pushed outwards. Indeed, transient red-shifted lines in the (5 Pictoris disks 
have been attributed to falling evaporating bodies (FEBs; Beust et al. 1990, 1996, 1998). 
While these events occur only for bodies that pass within 0.4 AU of the star, sublimation of 
comets and asteroids can occur as far out as 4 AU (Flammer et al. 1998). Since the dust/gas 
ratio in comets is close to 1 (Greenberg 1998), large amounts of evaporated fine dust par- 
ticles might exist in the region of 49 Cet's inner disk. Given the wide-ranging grain sizes 
in optically thin circumstellar disks, from submicron-sized grains in 49 Cet and HD 141569 
(Marsh et al. 2002) to the 50/im grains in part of the disk around HR 4796A (Wahhaj et al. 
2005), we conclude that multiple processes may control the evolution of dust grains both in 
size and radial location. 

In summary our analysis of 49 Cet, (3 Pic and HR 4796A demonstrates that their disks 
have more in common than just a large fractional luminosity. Morphological resemblances 
that set these objects apart from other "Vega-like" A stars, including a Psa and Vega itself 
may be due to their youth, 8-20 Myr (see discussion in Jura et al. 1993). They may in fact be 
in a transitional stage of evolution between that of viscous accretion disks around Herbig Ae 
stars (e.g., Mannings & Sargent 1997; Mannings, Koerner & Sargent 1997) and truly "Vega- 
like" stars with ages of order 100 Myr or greater. As such, they are interesting candidates 
for constraints on the survival times of circumstellar gas. Marginal detection of molecular 
gas (~ 10~ 3 M ) has been reported around 49 Cet (Thi et al. 2001) based on observations 
of the pure rotational transitions of H 2 . If circumstellar exists, it will be an important factor 
in the evolution of dust grains, since this amount of gas will shorten the natural lifetime even 
for grains as large as a ~ 50 fim (Klahr & Lin 2001). A marginal detection of CO(2^1) 
emission at 1mm was also reported by Zuckerman et al. (1995). However, observations of 
CO (3— >2) at 850 /xm by Coulson et al. (1998) suggest a very low upper limit on the gas mass 
in the disk < 6 x 10~ 8 M Q . CO is easily frozen out onto grains in the colder parts of the 
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disk (Aikawa et al. 1997). 

A better understanding of both the abundance of molecular gas and the morphology 
of cold dust in the outer disk await study at longer wavelengths with higher resolution and 
sensitivity, such as will be available to the next generation of mm-wave arrays. 

Observations presented here were obtained at the W.M. Keck Observatory (WMKO), 
which is operated as a scientific partnership among the California Institute of Technology 
the University of California and the National Aeronautics and Space Administration. The 
Observatory was made possible by the generous financial support of the W.M. Keck Foun- 
dation. A great debt is due, also, to Robert Goodrich and the WMKO summit staff for 
their many hours of assistance in adapting MIRLIN to the Keck II visitor instrument port. 
The authors wish also to recognize and acknowledge the very significant cultural role and 
reverence that the summit of Mauna Kea has always had within the indigenous Hawaiian 
community. We are most fortunate to have the opportunity to conduct observations from 
this mountain. 
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Table 1. Flux Densities for 49 Cet 
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Jayawardhana et al. 2001 


12.5 


1.2 


0.2 


0.026 


0.143 


This work 


17.9 


2.0 


0.186 


0.025 


0.069 




12.0 


6.5 


0.33 


0.066 


0.154 


IRAS FSC 


25.0 


11.0 


0.38 


0.076 


0.036 


» 


60.0 


40.0 


2.0 


0.4 


0.006 


5) 


100.0 


37.0 


1.91 


0.38 


0.002 


5) 


150.0 


40.0 


0.75 


0.50 


~ 


ISO 


170.0 


50.0 


1.1 


0.50 


~ 


5) 


800.0 




<0.036 (3a) 




~ 


Zuckerman & Becklin 1993 
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0.0127 
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Bockelee-Morvan(IRAM) 1994 
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Arcseconds 




Fig. 1.— KECK/MIRLIN contour maps of 49 Cet at A = 12.5 and 17.9 /jm. North is "up". 
Contour images of the PSF star, a Boo are displayed in panels a & b as insets with contour 
levels of 10%, starting at the 60% level. The FWHM of the PSFs are 0".41 and 0".48 at 
12.5 and 17.9 fim respectively, (a) Emission at 12.5 fim. Lowest contour level is at the 2a 
level (41 mJy/arcsec 2 ). Higher contour levels are at 2a x (io - 068 ™) for the n th contour, (b) 
Emission at 17.9 /im contoured as in a. Lowest contour is at the 2a level (26 mJy/arcsec 2 ). 
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Astrocentric Distance (AU) Astrocentric Distance (AU) 

Fig. 2. — Flux density strip cuts of width equal to FWHMs of PSFs taken along the long 
axes of emission ( PA = 125°) at 12.5 and 17.9//m seen in Fig. 1. Solid lines represent the 
49 Cet emission profiles and dot-dashed lines the profile for PSF star a Boo. Dashed lines 
are from best-fit models of emission described in Section 3 of the text. 



-16- 



Arcseconds 

-10 1-10 1 



100 



5 ° 



-100 



-50 50 -50 50 

AU 

Fig. 3. — Simulated images at 12.5 /xm and 17.9 /im from the best-fit model. Contour levels 
are the same as in Fig. 1. For this model, = 57°, 9 = 125°, R in = 35 AU, AR = 10 AU, 
a = 0.1 /im, and a = 1.3 x 10~ 3 for the inner disk, and = 60°, 9 = 125°, R in = 45 AU, 
AR = 900 AU, a = 23 /im, and a = 5.4 x 10 -4 for the outer disk. Since only very weak 
constraints were found for the 7 and (3 parameters they were set to the nominal values of 
and 1 respectively. 
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Fig. 4. — Flux densities of 49 Cet from Table. 1 (circles) and UBVJHKL & M photometry 
from Sylvester et al. 1996 (triangles). The solid line represents the complete SED for the best- 
fit model shown in Fig. [3j The dashed line displays the contribution from the photosphere 
derived from a 1993 Kurucz model fit to the short-wavelength photometry. The dash-dot 
line is the contribution from the inner disk suggested by the Keck MIR images. The dotted 
line is contribution from the outer disk. 
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Fig. 5. — Probability distributions for selected inner-disk model parameters. Most probable 
values correspond to peaks in the probability distributions. The uncertainties are estimated 
as the 66% confidence intervals (where our distributions are binned finely enough) or the 
shortest range of parameter values that encompasses 66% of the total probability. The most 
probable values and uncertainties for these parameters are detailed in the text. 
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Fig. 6. — Probability distributions for selected outer-disk model parameters. 



